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I.  INTRODUCTION 


In  the  past  40  years  considerable  effort  has  been  directed  toward  analyzing 
the  physical  and  chemical  processes  involved  in  solid  propellant  combustion. 

Part  of  this  effort  has  been  concerned  with  the  development  of  idealized  models 
of  the  combustion  process  as  a  whole  in  order  to  determine  how  the  magnitude, 
speed,  and  site  of  the  chemical  energy  release  controls  the  burning  rate.  With 
few  exceptions  a  common  thread  running  through  these  treatments  is  the  use  of 
the  experimental  burning  rate  or  its  derivatives  (i.e.,  the  pressure  and  tempera¬ 
ture  sensitivities)  as  the  authenticating  test  of  the  model.  Because  of  its 
long  chain  of  continuity  and  its  comparative  dominance  of  the  field  we  shall  term 
this  validation  scheme  the  Classical  Approach.  In  this  paper  we  shall  reexamine 
the  utility  and  limitations  of  this  approach  in  unfolding  solid  propellant  burning 
mechanisms . 

For  the  sake  of  clarity  our  usage  of  a  few  key  terms  should  be  explained.  We 
consider  that  a  solid  propellant  combustion  model  possesses  two  essential  features: 
an  idealization  of  the  burning  process,  i.e.,  a  set  of  differential  equations  and 
boundary  conditions  presumed  to  express  formally  the  combustion  mechanism,  and 
an  algorithm  or  computational  scheme  for  solving  these  equations.  The  word 
'’mechanism"  is  used  to  describe  the  collection  of  fundamental  physical  and  chemical 
processes  proposed  as  responsible  for  the  self-sustained  burning. 

The  most  commonly  assumed  mechanism,  and  the  one  to  which  we  shall  confine 
our  analysis  here,  consists  of  an  exothermic  surface  pyrolysis  reaction  which 
provides  the  reactants  for  a  single  exothermic  gas  phase  reaction.  One  formal 
idealization  corresponding  to  this  mechanism,  developed  in  detail  in  Reference 
1  involves  such  assumptions  as  Arrhenius  reaction  rates,  constant  specific  heats 
and  thermal  conductivity,  and  unity  Lewis  number.  These  assumptions  are 
typical  of  models  addressing  the  above  mechanism.  The  source  of  greatest  variety 
among  the  various  models  is  the  computational  algorithm, which  invariably  involves 
further  simplification  of  the  idealization. 

The  Classical  Approach  to  propellant  modeling  is  most  often  implemented  by 
developing  a  model  and  then  varying  the  kinetics  parameters  until  the  model  burning 
rates  closely  match  the  experimental  burning  rates.  Apart  from  limitations  in  the 
Classical  Approach  itself  (to  be  discussed)  this  particular  method  of  implementation 
raises  two  objections.  The  first  is  that  our  perception  of  the  phenomenology 
arising  from  the  assumed  mechanism  depends  on  the  fidelity  of  the  model  algorithms, 
and  these,  for  the  most  part,  have  not  been  tested.  Thus  it  is  possible  (see  the 
NC2  case^)  for  an  algorithm  to  be  particularly  inaccurate  for  the  set  of  parameters 
deduced  from  burning  rate  data  using  this  very  algorithm.  The  question  of  algorithm 
reliability  is  treated  in  a  companion  paper. 2  A  second,  more  subtle  objection 
to  this  method  arises  from  the  lack  of  uniqueness  associated  with  the  determined 
parameter  set.  The  functional  dependence  of  the  experimental  burning  rate  on 

^ Miller,  M.S.  ,  "In  Search  of  an  Idealized  Model  of  Homogeneous  Solid  Vropellant 
Combustion  >  "  Combustion  and  Flame,  Vol.  46,  pp.  51-73 ,  1982 . 

o 

Miller,  M.S.  and  Coffee,  T.P.,  "On  the  Numerical  Accuracy  of  Solid  Propellant 
Combustion  Models  ,"  accepted  for  publication  in  Combustion  and  Flam. 
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pressure  is  usually  quite  simple,  requiring  only  two  or  three  parameters  for  a 
close  fit  using  empirical  functions  such  as  r  =  a+bP  or  r  =  a+bPn.  Since 
typically  some  half  dozen  model  parameters  are  "floated",  it  is  an  open  question 
whether  a  successful  fit  is  due  to  the  correctness  of  the  assumed  mechanism  or 
simply  due  to  the  functional  flexibility  afforded  by  the  large  number  of  dispos¬ 
able  parameters.  In  view  of  the  fact  that  the  mechanism  being  considered  is 
known  (even  intended)  to  be  a  crude  representation  of  the  actual  phenomena  (in 
order  to  isolate  the  dominant  factors),  we  believe  that  an  exact  match  of  model 
to  experiment  is  simply  not  a  meaningful  objective. 

Recognizing  that  the  rationale  for  propellant  modeling  is  essentially 
heuristic,  one  might  be  better  advised  to  explore  thoroughly  the  range  of  burning 
rate  behavior  implicit  to  the  basic  idealization.  With  such  a  broad  understand¬ 
ing  of  the  model  phenomenology  one  could  then  attempt  to  identify  similar  trends 
in  real  propellant  combustion.  Such  a  strategy  is  clearly  incapable  of  validating 
the  proposed  mechanism;  however,  we  contend  that  this  limitation  is  intrinsic 
to  the  Classical  Approach  itself.  Further  progress  in  the  determination  of  the 
physical  mechanisms,  in  our  opinion,  can  only  result  from  more  detailed  experi¬ 
mental  constraints  on  the  structure  of  the  combustion  wave.  The  Classical 
Approach  is  no  more  and  no  less  than  a  preliminary  test  of  the  internal  consist¬ 
ency  and  compatability  of  potentially  useful  concepts. 

The  central  task  of  this  paper  then  is  to  explore  the  range  of  behavior 
of  the  burning  rate  expected  of  the  stated  idealization  using  a  numerically 
accurate  algorithm.  This  undertaking  is  complicated  by  the  possibility  that  the 
burning  phenomenology  may  depend  on  the  values  assigned  to  the  kinetics  parameters 
To  overcome  this  obstacle  we  have  put  together  a  number  of  different  data  sets 
which  cover  a  wide  range  of  parameter  combinations.  We  believe  that  the  range 
is  sufficiently  great  to  reveal  thoroughly  the  behavior  patterns  inherent  to  the 
mechanism.  One  can  think  of  these  data  sets  coupled  with  the  idealization  as 
model  propellants.  Following  the  Classical  Approach  we  shall  compare  the  way 
these  model  propellants  burn  with  real  propellants  and  attempt  to  draw  parallels. 

II.  THE  MODEL  PROPELLANTS 

Five  data  sets  will  be  developed  in  this  section  by  several  different  methods 
RDX1  and  NCI,  suggestive  of  nitramine  and  nitrate  ester  propellants,  are  obtained 
by  making  plausible  assumptions  and  estimates  based  largely  on  experiments  not 
involving  steady-state  burning.  NC2  is  based  on  the  kinetics  inferred  by  the 
method  of  fitting  a  model  to  double-base  rate  data.  The  self-consistency  of  this 
procedure  is  examined  elsewhere we  use  these  parameters  here  simply  to  explore 
as  wide  a  range  of  kinetics  as  possible.  NC3  is  assembled  almost  entirely  from 
Zenin 1 s  analysis^  of  embedded  thermocouple  data  taken  under  actual  burning 
conditions.  T1  was  put  together  as  a  test  case  to  exploit  the  availability  of  an 
exact  analytic  solution  found  for  the  heat  feedback  under  certain  conditions  • ^ 
This?test  case  also  served  as  a  partial  check  on  the  numerical  code  used  in  this 
work2.  The  parameter  values  determined  for  each  of  the  model  propellants  are 
summarized  in  Table  1. 


^ Zenin,  A. A.,  "Formal  Kinetic  Characteristics  of  the  Reactions  Accompanying  the 
Burning  of  a  Powder,  "  Fizika  Gorenii^a  i  Vzryva,  Vol.  2,  pp.  28-32,  1966. 
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TABLE  I. 


RDX1 

NCI 

NC2 

NC3 

T1 

2 

Mo(gm/cm  sec) 

1.035E12 

8.0E10 

3456 

460 

1.0E12 

Es (kcal/mole) 

47.8 

42 

10 

8 

45.0 

Qs(cal/gm) 

160 

75 

53 

75 

75 

ag 

2.5E16 

(sec-1) 

2.5E9 

(cc/gm  sec) 

1.44E12 

(sec-1) 

1.2E5 

(sec-1) 

100T 

(sec-1) 

E„ (kcal/mole) 

46.2 

•  15 

54 

5 

0 

QG(cal/gm) 

460 

400 

422 

350 

500 

WB(gm/mole) 

222 

40 

40 

40 

50 

V 

1 

2 

1 

1 

1 

-0 

mB 

.48 

1 

1 

1 

1 

N2 

6 

1.33 

1.33 

1.33 

1 

3 

Ps(gm/cm  ) 

1.6 

1.6 

1.55 

1.6 

1.6 

X  =  2.0E  - 

4  cal/cm 

sec°K  C 

P 

=  .35  cal/gm°K 
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A.  NCI 


The  rate  of  slow  thermal  decomposition  of  nitrocellulose  (NC)  and  NC-based 
propellants  has  been  measured  by  a  number  of  methods  including  weight  loss4, 
infrared  signature4,  evolution  of  gas  in  solvents5,  and  pulsed  calorimetry. 
Expressed  in  Arrhenius  form  the  frequency  factors  range  from  lO^-8  to  10^9  sec” 
and  the  activation  energies  from  42  kcal  to  46  kcal.  Thompson  and  Suh7  used 
an  ignition  delay  technique  and  obtained  a  frequency  factor  of  10^7sec-^  with 
an  activation  energy  of  40  kcal/mole.  As  representative  of  these  measurements 
we  take  the  frequency  factor  (A  )  to  be  7.5  X  lO^7  sec“l  and  the  activation 
energy  (Es)  to  be  42  kcal/mole. s  This  rate  is  used  for  the  controlling  solid 
reaction  and  is  assumed  to  occur  only  at  the  surface  of  the  burning  propellant. 
Since  this  activation  energy  is  suggestive  of  the  endothermic  scission  of  NO2, 
we  shall  assume  that  subsequent  exothermic  polymer  breakdown  reactions  occur 
quickly  on  the  surface  leading  to  a  net  exothermic  heat  release  Q  .  The  embedded 
thermocouple  measurements  of  Kubota,  et.al.8  and  Zenin9  indicate  a  value  of  Q 
dependent  on  the  burning  rate  in  the  range  60-130  cal/gm.  We  assume  the  constant 
value  75  cal/gm.  In  order  to  obtain  a  surface  pyrolysis  law  in  the  customary 
form 

M  =  Mq  exp(-Es/RTs)  (1) 

Mq  is  constructed  as  the  product  of  the  frequency  factor,  the  number  of  molecules 
per  unit  surface  area,  and  the  mass  per  molecule  (W.) . 

A 


M 

o 


The  propellant  mass  density  is 


and  N 

o 


is 


Avogadro’s  number. 


(2) 


4 

Phillips,  R.W.,  Or  lick,  C.A.  and  Steinberger,  R. ,  "The  Kinetics  of  the  Thermal 
Decomposition  of  Nitrocellulose,"  J_.  Ffo£s_.  CherrUj  Vol.  59,  p.  1034-1039,  1955. 

** Smith ,  R.D. ,  "Pyrolysis  of  Dissolved  Nitrocellulose,"  Nature,  Vol.  170,  p.  844- 
845,  1952. 

Q 

Aleksandrov,  V.  V. ,  Bufetov,  N.S.,  Pastukhova,  T.V.  and  Tukhtaev,  R.K.,  "Use 
of  Pulsed  Calorimetry  for  Investigating  the  Kinetics  of  Reactions  in  Condensed 
Mediai  Fizika  Goreniya  i  Vzrtjva,  Vol.  9,  pp.  75-83,  1973. 

7 

Thompson,  C.L.  and  Suh,  N.P.,  "The  Interaction  of  Thermal  Radiation  and  M-2 
Double-Base  Solid  Propellant, "  Comb_.  Sci_.  Tech. ,  Vol.  2,  pp.  59-66,  1970. 

g 

Kubota,  N.,  Ohlemiller,  T.J.,  Caveny,  L.H.  and  Summerfield,  M. ,  "The  Mechanism 
of  Super-Rate  Burning  of  Catalyzed  Double-Base  Propellants, "  AMS  Report  No.  1087, 
Princeton  University,  March  1973. 

g 

Zenin,  A. A.,  "Structure  of  Temperature  Distribution  in  Steady-State  Burning  of 
a  Ballistic  Powder j "  Fizika  Goreniya  i  Vzryva,  Vol.  2  pp.  67-76,  1966. 
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There  is  evidence  that  aldehyde-NC>2  reactions  may  play  a  significant  role 
in  the  fizz  zone  energy  release.  In  the  NCI  data  set  we  assume  they  are  dominant. 
Although  the  rates  of  these  reactions  are  not  precisely  known,  we  take  a  frequency 
9  cc 

factor  of  2.5  X  10  -  and  an  activation  energy  of  15  kcal/mole  to  be 

j  s  0  c 

representative.  The  values  of  other  necessary  parameters  in  Table  1  are  reason¬ 
ably  consistent  with  this  assumed  reaction  identity  although  the  nascent  reactant 
mass  fraction  (mg-®)  is  assigned  arbitrarily.  Values  for  the  thermal  conductivity 
and  specific  heat  are  assumed  to  be  the  same  for  all  of  the  model  propellants  and 
are  estimated  from  Zenin.® 


B.  RDX1 


The  most  comprehensive  attempt  to  model  RDX  monopropellant  combustion  is 
due  to  BenReuven,  et  This  model  assumes  that  the  surface  regresses 

primarily  by  RDX  evaporation  from  a  melt  layer.  Two  gas -phase  reactions  are  then 
considered;  unimolecular  decomposition  of  RDX  vapor  and  formaldehyde-NC^  reactions. 
The  latter  are  found  to  be  responsible  for  only  a  small  percentage  (^15%  at  40atm) 
of  the  heat  feedback.  Although  exothermic  decomposition  of  RDX  is  allowed  in  the 
melt  layer,  it  appears  that  the  net  effect  of  these  processes  is  endothermic. 

Very  recent  embedded  thermocouple  results^  for  RDX  composite  propellants 
(^80%  RDX) ,  however,  indicate  an  even  higher  exothermicity  to  the  surface  re¬ 
action  than  for  NC-based  propellant. 

Using  these  observations  to  simplify  the  combustion  mechanism  to  only  two 
reactions,  we  have  developed  the  following  conjectural  picture  of  RDX  combustion. 
When  an  RDX  molecule  at  or  very  near  the  surface  decomposes,  the  chemical  energy 
which  becomes  available  may  be  used  in  either  of  two  ways.  One  or  more  unreacted 


40Lengelle,  G. ,  Duterque,  J . ,  Verdier,  C. ,  Bizot,  A.  and  Trubert,  J. ,  "Combustion 
Mechanisms  of  Double  Base  Solid  Propellants ,"  Seventeenth  Symposium  (Inters 
national)  on  Combustion,  The  Combustion  Institute,  pp.  1443-1451 ,  1978. 

11  Pollard,  F.H.  and  Wyatt,  "Reactions  Between  Formaldehyde  and  nitrogen 

Dioxide  I.,"  Trans_.  Faradmj_  Soc. ,  Vol.  45,  pp.  760-767,  1949. 

12 BenReuven,  M. ,  Caveny.  L.J.,  Vichnevetsky ,  R.J.,  Summer field,  M.,  "Flame  Zone 
and  Sub-Surface  Reaction  Model  for  Deflagrating  RDX,"  Sixteenth  Symposium 
(International)  on  Combustion,  The  Combustion  Institute,  pp.  1223-1233,  1977. 

BenReuven,  M. ,  and  Caveny,  L.H.,  "Nitramine  Flame  Chemistry  and  Deflagration 
Interpreted  in  Terms  of  a  Flame  Model,  AIAA  Paper  79-1133,  AIAA/SAE/ASME 
15th  Joint  Propulsion  Conference  (Las  Vegas,  NV),  June  18-20 ,  1979. 

14 BenReuven,  M. ,  and  Caveny,  L.H.,  "Nitramine  Monopellant  Deflagration  and 
General  Nonsteady  Reacting  Rocket  Chamber  Flows,  '  MAE  Report  No.  1455 , 
Princeton,  University,  January  1980. 

15 Kubota,  N. ,  "Combustion  Mechanisms  of  Nitramine  Composite  Propellants, " 

Eighteenth  Symposium  (International)  on  Combustion,  The  Combustion  Institute, 
pp.  187-194,  1981. 
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molecules  may  use  this  energy  to  escape  the  surface  bonding  forces  into  the  gas 
phase,  or  the  energy  may  be  absorbed  into  the  surface  as  heat.  In  the  absence 
of  more  detailed  knowledge  we  assume  a  simple  statistical  partitioning  of  the 
energy  among  these  possible  modes.  If  we  assume  with  BenReuven,  et  al12  that 
RDX  decomposes  into  the  following  intermediates 

RDX  |  N2  +  N20  +  N02  +  3CH20 


then  the  enthalpy  change  (solid  RDX  to  gaseous  products)  at  298°K  is  317  cal/gm 
of  RDX.  (Heat  of  formation  for  RDX  is  from  Dobratz16  and  that  for  the  products 
from  the  JANNAF  tables17)  ■  Assuming  that  half  of  the  317  cal/gm  is  used  to  eject 
unreacted  RDX  molecules  and  half  goes  to  surface  heating,  we  obtain  Q  =  160  cal/gm 
Using  the  heat  of  sublimation  measured  by  Cundall1^  (145  cal/gm)  as  aSmeasure  of 
the  surface  bonding  energy  of  RDX  molecules,  we  compute  that  1.09  gm  (=317/2/145) 
of  RDX  is  ejected  unreacted  into  the  gas  phase  per  gram  of  RDX  that  undergoes 
surface  pyrolysis.  Using  Eq.  (2)  and  noting  that  the  surface  recedes  as  a 
result  of  RDX  loss  by  pyrolysis  and  ejection. 


M 

o 


(1+1.09) 


1/3 


A 

s 


gm 

2 

cm  sec 


(3) 


The  condensed  phase  kinetics  (specifically,  the  liquid  phase)  for  RDX  are  taken 
from  the  review  by  Schroeder.19  To  be  consistent  with  this  use  of  liquid  phase 
kinetics,  an  explicit  treatment  of  the  melt  layer  in  our  energy  partitioning 
calculation  might  be  desirable;  however  the  above  estimates  are  adequate  to  the 
purpose  of  this  paper. 


After  surface  decomposition  by  what  might  be  termed  ’’dispersive  pyrolysis", 
it  is  assumed  that  only  the  vapor  phase  decomposition  of  the  unreacted  RDX 
results  in  significant  heat  feedback  to  the  surface.  This  rate  is  also  taken 
from  Schroeder.1^  The  gas-phase  reaction  heat  is  obtained  by  summing  the 
heat  of  sublimation  and  the  reaction  enthalpy  of  solid  RDX  to  gaseous  products. 


16  // 
Dobratz,  B.M.,  "Properties  of  Chemical  Explosives  and  Explosive  Simulants t 

University  of  California  (Livermore)  Report  No.  UCRL-51319,  Rev.  1 ,  July, 

1974. 

1 7 

Stull,  D.R.  and  Prophet,  H.  JANAF  Thermochemical  Tables ,  2nd  Edition, 
NSRDS-NBS-37,  June  1971. 

1 8 

Cundall,  R.B.,  Palmer,  T.F.,  Wood,  C.E.C.,  "Vapor  Pressure  Measurements  on 
Some  Orqanic  Eiqh  Explosives ,  J.  Chem.  Soc.,  Faraday  Trans.,  Vol.  I, 
pp.  1339-1345,  1978. 

19 

Schroeder,  M.A.,  "Critical  Analysis  of  Nitramine  Decomposition  Data:  Activa¬ 
tion  Energies  and  Frequency  Factors  for  HMX  and  RDX  Decomposition,  Proceedings 
of  17th  JANNAF  Combustion  Meeting,  CPIA  Pub.  No.  329,  Vol.  II,  pp.  493-508, 
November  1980. 
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C.  NC  2 


20 

Beckstead  has  recently  applied  the  Classical  Approach  to  double  base 
propellants  of  varying  heats  of  explosion (HEX) .  In  order  to  achieve  a  fit  of 
the  standard  BDP  monopropellant  model ^  to  the  burning  rate  data  at  a  given  HEX, 
the  model  "flame11  temperature  Tf  (actually  the  dark  zone  temperature)  was  allow¬ 
ed  to  vary  with  pressure.  Thus  at  a  given  pressure  the  usual  BDP  formalism 
applies.  Values  of  Es  and  E q  were  determined  as  10  kcal/mole  and  54  kcal/mole, 
respectively.  Not  all  of  the  other  fitted  parameters  were  reported  but  we  have 
used  what  information  was  given  to  construct  a  consistent  data  set  for  HEX=1000 
cal/gm  at  136  atm,  the  highest  pressure  considered.  In  order  to  do  this  we 
first  assumed  the  values  of  m^-0,  W^,  and  shown  in  Table  1. 


D.  NC3 

3  9  22  23 

Zenin  9  9  9  has  reported  extensive  measurements  and  analysis  of 

temperature  distributions  in  burning  double-base  propellants  obtained  by  the 
embedded  thermocouple  technique.  Analysis  of  the  spatial  heat  release  in  the 
condensed  phase  suggests  a  reactive  zone  of  at  least  10  microns  thickness. 
However,  in  order  to  use  the  Eq.  (1)  form  of  the  pyrolysis  law,  we  have  con¬ 
structed  an  Arrhenius  plot  of  M  vs.  1/TS  and  determined  values  of  M0  and  Es 
consistent  with  Zenin' s  measurements . 0,23  As  in  NCI,  we  have  taken  Qs  to  be 
75  cal/gm. 

The  fizz  zone  frequency  factor  was  deduced  as  follows.  Values  of  1,40, 
and  30  were  assigned  to  mg'"0,  Wg,WQ,  respectively.  From  Reference  9  a  value 
of  Qg  =  350  cal/gm  was  selected.  Zenin3  determined  the  fizz  zone  activation 
energy  to  be  5  kcal/mole  and  the  reaction  order  to  be  one.  We  then  solve  the 
expression^ 

q  =  Cp(T£-T)(^)AG  exp(-EG/RT)  (4) 

for  Aq  at  the  peak  in  the  fizz  zone  energy  release  profile  shown  in  Fig.  1  of 
Reference  3. 


20 Beckstead,  M.W.,  " Model  for  Double-Base  Propellant  Combustion, '  AIM  Journal , 
Vol.  18,  pp.  980-985,  1980. 

21 Beckstead,  M.  W. ,  Derr,  R. L.  and  Price,  C.  E. ,  "The  Combustion  of  Solid  Mono¬ 
propellants  and  Composite  Propellants Thirteenth  Symposium  (International) 
on  Combustion,  The  Combustion  Institute,  pp.  1047-1056,  1971. 

22 Zenin,  A. A.,  "Burning  of  Nitroglycerine  Powder  in  Vaccuum  and  at  Subatmo spheric 
Pressures,"  Fizika  Gorenipya  i  Vzryva,  Vol.  2,  pp.  74-78,  1966. 

23 Zenin,  A. A.  and  Nefedova,  0.1. .  "Burning  of  Ballistic  Powder  Over  a  Broad 
Range  of  Initial  Temperatures/'  Fizika  Goreniya  i  Vzryva,  Vol.  5,  pp.  45-53, 
1967. 
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E.  T1 


The  test  case  T1  is  assigned  values  for  the  pyrolysis  law  in  the  same 
general  range  as  NCI  and  RDX1.  The  temperature  dependence  of  the  frequency 
factor  for  the  gas  reaction  is  chosen  so  that  the  reaction  rate  depends  only 
on  the  local  mass  fraction  and  not  on  local  temperature.  This  enables  use  of 
an  exact  solution  for  the  heat  feedback  (Eq. (31)  of  Ref.  1).  The  range  of 
values  for  unimolecular  reactions  is  large^4,25  and  the  chosen  rate  lies  well 
within  this  range.  The  remaining  parameter  values  are  chosen  to  be  similar 
to  the  other  data  sets. 


III.  REAL  VS.  MODEL  PROPELLANT  COMBUSTION 


In  this  section  we  shall  examine  the  combustion  patterns  peculiar  to  the 
chosen  idealization  and  determine  the  extent  to  which  these  patterns  parallel 
real  propellant  combustion.  Though  many  calculations  of  burning  rate  based  on 
this  idealization  have  been  published  we  contend  that  these  patterns  have  not 
been  reliably  identified  due  to  the  limitations  of  using  a  single  data  set  or 
an  inaccurate  algorithm  or  both.  Here,  with  the  exception  of  the  NC2  and  T1 
data  sets,  the  model  propellant  burning  rates  were  calculated  by  accurate 
numerical  integrations."  The  asymptotic  theory  of  Williams2  modified  by  the 
adiabatic  deflagration  limits  of  Buckmaster,  et  al22  was  used  for  NC2  after 
validation2  by  numerical  integration  at  136atm  for  T0  in  the  range  -123°C  to 
+77°C .  The  heat  feedback  for  the  T1  case  could  be  obtained  analytically  (using 
Eq.  31  of  Ref.  1)  so  that  the  burning  rate  determination  is  simply  algebraic 
(although  non-analytic) .  The  variation  of  these  model  burning  rates  with 
pressure  and  initial  temperature  is  shown  in  Figs.  1-10. 


The  general  pressure  dependence  of  the  model  burning  rates  is  found  to  be 
well  approximated  by  the  expression  r=a+bP  where  n^is  slightly  less  than  half 
the  gas  reaction  order  (v) .  As  discussed  elsewhere  the  idealization  considered 
leads  to  a  constant  rate  at  sufficiently  low  pressure  where  the  gas-phase  heat- 
feedback  is  negligible  compared  with  the  solid-phase  heat  release.  With  a 
laminar  gas  flame  the  pressure  index  is  exactly  v/2  so  that  the  solid  regression 
rate  has  much  the  same  character  over  at  least  some  pressure  range.  In  general 
the  higher  the  gas  reaction  activation  energy  the  closer  n  is  to  v/2.  The  high 


24 Adams ,  G.F.,  "An  Analysis  of  the  Pressure  Dependence  of  Nitrate  Ester  Thermal 
Decomposition,  "  Ballistic  Research  Laboratory  Technical  Report  ARBRL-TR-02106 , 
September,  1978 . 

25 Adams,  G.F.,  "A  Priori  Estimation  of  Rate  Constants  for  Unimolecular  Decomposi¬ 
tion  Reactions /'  Ballistic  Research  Laboratory  Technical  Report  ARBRL-TR-02143, 
February,  1979. 

26 Williams,  F.A.,  "Quasi-Steady  Gas-Flame  Theory  in  Unsteady  Burning  of  a 
Homogeneous  Solid  Propellant AIAA  Journal,  Vol.  11,  pp.  1328-1330,  1973. 

27 

Buckmaster,  J.D.,  Kapila,  A.K.,  and  Ludford,  G.S.S.,  "Linear  Condensate 
Deflacrnation  for  Large  Activation  Energy,"  Acta  A^tronautica,  Vol.  3, 
pp.  593-614,  1976. 
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Figure  1.  Burning  rate  vs.  pressure  computed  for  the  RDX1 
data  set  at  various  initial  temperatures. 
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Fig.  2.  Burning  rate  ys,  pressure  computed  for  the  NCI 
data  set  at  various  initial  temperatures . 
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Fig.  3.  Burning  rate  ys.  pressure  computed  for  the  NC2  data  set 
at  various  initial  temperatures. 
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Fig.  4.  Burning  rate  vs.  pressure  computed  for  the  NC3 
data  set  at  various  initial  temperatures. 
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Fig.  5.  Burning  rate  vs.  pressure  computed  for  the  T1 
data  set  at  various  initial  temperatures. 
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Fig.  6.  Burning  rate  vs.  initial  temperature  computed  for  the 
RDX1  data  set  at  various  pressures. 
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Fig.  7.  Burning  rate  vs.  initial  temperature  computed  for 
the  NCI  data  set  at  various  pressures. 
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Fig.  8.  Burning  rate  vs.  initial  temperature  computed  for  the 
NC2  data  set  at  various  pressures. 
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Fig.  9.  Burning  rate  vs.  initial  temperature  computed  for 
the  NC3  data  set  at  various  pressures. 
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Fig.  10.  Burning  rate  vs.  initial  temperature  computed  for  the 
T1  data  set  at  various  pressures. 
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pressure  deflagration  limit  reached  when  essentially  all  of  the  gas  heat  release 
is  returned  to  the  surface  has  no  practical  effect  on  the  model  propellant  burn 
rates.  Its  influence  can  be  detected  however  in  the  slight  pressure  index 
decrease  between  100  and  1000  atm  for  NC3  at  the  lowest  T0.  For  some  other 
combination  of  parameters  its  effect  may  be  greater.  Note  that  for  several  of 
the  model  propellants  no  departure  of  the  rate  from  P  dependence  is  observed 
over  the  chosen  1-1000  atm  range  even  though  the  value  of  n  does  decrease 
slightly  with  increasing  T  . 

The  initial  temperature  dependence  of  the  model  burning  rates  is  displayed 
as  Zn  r  vs.  T0  to  facilitate  comparison  with  the  real  propellant  burning  rates 
to  be  discussed.  This  dependency  is  usually  discussed,  however,  in  terms  of 
the  temperature  sensitivity.  Op,  defined  as  the  slope  of  the  Zn  r  vs.  T0  curve. 

The  solid  curve  in  Fig.  11  summarizes  the  general  shape  of  Op(T0)  found  for  the 
model  propellants.  At  a  higher  pressure  the  curve  shifts  to  the  right  as  in¬ 
dicated  by  the  dashed  line.  (op(T0)  is  plotted  explicitly  for  each  data  set 
in  Ref.  2).  Any  given  model  propellant  will  exhibit  only  part  of  this  charac¬ 
teristic  curve  over  the  practical  range  of  T0  and  P.  The  shape  of  crp(T0)  in 
Fig.  11  arises  from  the  transition  between  a  laminar  flame  character  at  low  Tq 
to  a  rate  which  is  completely  controlled  by  the  heat  from  the  solid  reaction 
at  high  T0.  Movement  on  the  curve  through  the  transition  region  from  left  to 
right  indicates  a  decrease  in  the  contribution  to  the  surface  temperature  coming 
from  the  gas.  This  occurs  when  increases  in  Ts  (due  to  increases  on  Tq)  result 
in  a  lower  gas-phase  heat  feed-back.  This  happens  when  the  increase  in  mass 
flow  accompanying  the  Ts  increase  overcomes  the  increased  gas  reaction,  blowing 
the  flame  further  from  the  surface. 

In  order  to  compare  these  model  propellant  trends  to  real  propellants  we 
have  collected  representative  data  for  pure  nitramines  ("monopropellants")  and 
nitrate  ester  propellants  in  Figs.  12-19.  For  ready  comparison  these  data  are 
plotted  on  the  same  scales  used  for  the  model  propellants. 

As  with  the  model  propellants  the  real  propellant  burning  rates  can  be  well 
represented  by  the  form  r=a+bPn.  The  maximum  pressure  index  is  about  .8  -  .9 
compared  to  the  model  values  near  .5  or  1.  Also,  like  the  model  propellants, 
there  is  a  slight  decrease  in  the  value  of  the  pressure  index  at  any^given 
pressure  as  the  initial  temperature  is  raised.  Taking  20  atm  and  25°C  as  an 
arbitrary  reference,  the  real  burning  rates  are  closely  grouped  at  .3  -  .45 
cm/sec  whereas  the  model  propellants  span  the  wider  range  .1-8  cm/sec. 

At  20  atm  and  25°C  the  nitrate  ester  propellants  have  a  crp  of  about  .0056  - 
0059°C_1  and  the  value  for  HMX  is  about  . 0023°C"  .  The  model  propellant  Op  is 
in  the  range  .0023  -  .0076°C_1  at  the  same  temperature  and  pressure.  Over  Changes 
of  about  50°C  little  change  in  op  occurs  for  both  real  and  model  propellants. 

The  Zenin  data23  (Fig.  19)  cover  the  widest  range  of  T0  and  show  ap  increasing 
with  increasing  T0.  Other  nitrate  ester  data  (Fig.  17  and  18)  show  a  nearly 

constant  a  over  the  measured  range. 
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Fig.  11.  Patterns  of  burning  rate  control  evidenced  by  temperature 
sensitivity  calculations  for  idealized  propellants. 
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Fig.  12.  Reported  measurements^  ^  of  burning  rate  vs.  pressure 
for  RDX  monopropellant. 


27 


H  MX 

PRICE  ,  et  al.  ( 1979  ) 


—  1 
v% 

\ 

6 


'/ 

,'F 


°X 


.1  - — 


□  15  0°  C 
O  25°  C 


J _ I 


10 


P (atm) 


100 


Fig.  13. 


31 

Reported  measurements  of  burning  rate  vs.  pressure 
for  HMX  monopropellant. 
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Fig.  14.  Reported  measurements  of  burning  rate  vs.  pressure  for 
a  double-base  propellant  (coded  BC-EI-4  in  Ref.  32). 
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Fig.  15.  Reported  measurements^  of  burning  rate  vs.  pressure 
for  a  double-base  propellant  (N  powder). 
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Fig.  16.  Reported  measurements  ’  of  burning  rate  vs.  pressure 
for  a  double-base  propellant  (JP76). 
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Fig.  17.  Reported  measurements  of  burning  rate  vs.  initial 

temperature  for  a  double-base  propellant  (coded  BC-EI-4 
in  Ref.  32). 
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Fig.  18.  Reported  measurements  ’  of  burning  rate  vs.  initial 
temperature  for  a  double-base  propellant  (JP76). 
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Fig.  19.  Reported  measurements23  of  burning  rate  vs.  initial 
temperature  for  a  double-base  propellant  (N  powder). 
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This  comparison  of  real  and  model  propellant  phenomenology  can  be  used  to 
suggest  likely  weaknesses  in  the  idealization  under  consideration.  The 
discrepancy  in  maximum  pressure  index  might  arise  from  several  sources.  Thermo¬ 
couple  data^’P  have  indicated  that  both  the  solid  phase  heat  release  (Qs)  and 
the  dark  zone  temperature  (Tf)  increase  with  pressure.  Either  effect  could 
alter  the  effective  pressure  index  although  neither  would  affect  the  qualitative 
behavior  of  Op  in  Fig.  11.  The  assumption  of  constant  conductivity  and  specific 
heats  is  also  only  a  rough  approximation  and  could  influence  the  pressure 
index.  The  assumption  of  a  single  elementary  reaction  in  each  phase  is  only 
a  convenient  idealization.  The  presence  of  parallel  reactions  of  different 
orders  might  cause  the  pressure  index  to  increase  with  increasing  pressure. 

It  would  seem,  however,  that  little  would  be  gained  by  adding  any  of  these 
refinements  unless  they  first  could  be  characterized  through  experiments. 

IV.  CONCLUSIONS 

We  have  described  a  systematic  study  of  the  phenomenology  associated  with  a 

typical  idealization  of  solid  propellant  combustion.  Many  of  the  trends 

identified  for  the  "model  propellants"  were  shown  to  have  parallels  in  real 

propellant  combustion.  The  general  pressure  dependence  of  the  burning  rate 

for  both  model  and  real  propellants  is  similar  but  has  too  little  structure 

to  be  of  much  heuristic  value.  The  temperature  sensitivity  as  a  function  of 

initial  temperature  for  various  pressures  has  a  richer  structure,  the  details 

of  which  correlate  with  the  relative  control  of  the  burning  rate  by  solid  and 

gas  phase  processes.  Various  experimental  a  values  follow  different  features 

of  the  model  propellant  pattern,  but  in  mostPcases  data  over  a  wider  range  of 

T  are  needed  to  confirm  the  parallels, 
o 

The  Classical  Approach  to  propellant  modeling,  if  properly  applied,  can 
be  a  useful  method  of  exploring  the  cooperative  effects  of  proposed  combustion 
mechanisms.  An  important  limitation  to  its  use,  however,  derives  from  the 
fact  that  errors  in  the  mechanism  description  often  have  a  weak  effect  on  the 
burning  rate.  An  erroneously  high  reaction  rate  in  the  gas  phase,  for  example, 
tends  to  increase  the  heat  feedback  which  leads  to  a  higher  burning  rate.  But 
the  magnitude  of  the  burning  rate  increase  is  moderated  by  an  increased  con¬ 
vective  heat  loss.  The  net  effect  is  that  the  burning  rate  may  not  be  a  sensi¬ 
tive  function  of  the  mechanisms  on  which  it  depends.  For  this  reason  the  next 
step  in  modeling  propellant  combustion,  in  our  opinion,  will  have  to  be  based 
on  mechanisms  independently  validated  for  each  phase.  Such  mechanisms  need 
not  be  greatly  increased  in  sophistication  or  chemical  specificity  but  must 
be  based  on  characterizations  of  the  solid  and  gas  phase  processes  which  are 
consistent  with  experimental  measurements  of  the  detailed  combustion  wave 
structure. 
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A_  -  gas-phase  reaction  frequency  factor 
G 

-  condensed-phase  reaction  frequency  factor 

B  -  gas-phase  reactant  label 

C  -  specific  heat  for  solid  and  gas  phases 
P 

E_  -  activation  energy  for  gas-phase  reaction 
G 

E  -  activation  energy  for  solid-phase  reaction 
s 

M  -  mass  regression  rate  (mass  flux) 

M  -  constant  in  pyrolysis  surface  decomposition  mechanism 

m  “°,  m  +°  -  mass  fraction  of  B  evaluated  at  negative  and  positive  sides  of 
B  B  the  solid/gas  interface,  respectively. 

No  -  Avogadro's  number 

N0  -  number  of  moles  of  C  produced  per  mole  of  B  which  reacts 
n  -  pressure  exponent  of  burning  rate 
P  -  total  pressure 

Q  -  heat  of  reaction  per  unit  mass  for  solid  reaction  (positive  for  exothermic) 

Q  -  exothermic  gas-phase  reaction  heat  per  unit  mass  of  B 
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R  -  universal  gas  constant 
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T  T  T  -  initial,  surface,  and  flame  temperatures,  respectively 
o’  s’  f 

W.,WD,W„  -  molecular  weights  of  A,B,  and  C 
A  B  L 

W  -  average  molecular  weight  of  mixture  in  gas-phase 
X  -  heat  conductivity  in  gas-phase 
v  -  gas-phase  reaction  order 
Ps  -  mass  density  of  solid 

a  -  temperature  sensitivity  at  constant  pressure 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number^ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name  : 


Telephone  Number: 
Organization  Address: 


